The Role of Surface Plasmons in The Casimir Effect 



F. Intravaia/ C. Henkel/ and A. Lambrecht^ 

^ Institut fUr Physik, Universitdt Potsdam, Am Neuen Palais 10, 14469 Potsdam, Germany 
^Laboratoire Kastler-Brossel ENS, UPMC, CNRS, Case 74, Campus Jussieu, 75252 Paris Cedex 05, France 

In this paper we study the role of surface plasmon modes in the Casimir effect. The Casimir energy can be 
written as a sum over the modes of a real cavity and one may identify two sorts of modes, two evanescent surface 
plasmon modes and propagative modes. As one of the surface plasmon modes becomes propagative for some 
choice of parameters we adopt an adiabatic mode definition where we follow this mode into the propagative 
sector and count it together with the surface plasmon contribution, calling this contribution "plasmonic". We 
evaluate analytically the contribution of the plasmonic modes to the Casimir energy. Surprisingly we find that 
this becomes repulsive for intermediate and large mirror separations. The contribution of surface plasmons to the 
Casimir energy plays a fundamental role not only at short but also at large distances. This suggests possibilities 
to taylor the Casimir force via a manipulation of the surface plasmons properties. 

PACS numbers: 42.50.Pq Cavity quantum electrodynamics - 73.20.Mf Collective excitations 



I. INTRODUCTION 

The Casimir force is the archetypal mechanical conse- 
quence of vacuum fluctuations in the quantized electromag- 
netic field. In its simplest form, it gives rise to the attraction of 
two planar mirrors placed in empty space at zero temperature 
fl |. The corresponding interaction energy E takes a universal 
form for perfect reflectors, 

^ = ^^^^ = "4?^' 
where L is the distance between the mirrors, A their area, and 
h and c the reduced Planck constant and the speed of light. We 
abbreviate K = ISO/tt^ ^ 5.8052762. As usual in thermody- 
namics, a negative energy corresponds to a binding energy. 

The Casimir force was soon observed in different exper- 
iments which confirmed its existence El [3] |4l. In recent 
years, technological improvement allowed to reach a preci- 
sion in the percent range, which makes an accurate compar- 
ison to theoretical predictions possible and has prompted a 
series of refined calculations |5, 6|. Casimir's 1948 deriva- 
tion of Eq.Q is based on summing the zero-point energies 
^huj of the cavity eigenmodes, taking the difference for fi- 
nite and infinite separation, and removing the divergences by 
inserting a high-energy cutoff. He considered an ideal set- 
ting with perfectly reflecting mirrors in vacuum. Experiments 
are however performed with real reflectors, typically metallic 
mirrors which are good reflectors only at frequencies below 
the plasma frequency (cjp/27r) or alternatively at wavelengths 
much larger than Ap = 27r c/ujp. It has been known since 
a long time that this has a significant effect on the force, in 
particular at mirror distances of the order of Ap or smaller 
CI El 05 and precise investigations have been developed re- 
cently EiiaiiiiiiainiiTiiisiiisiiniii^ 

A system made from real material mirrors sustains electro- 
magnetic modes which strongly differ with respect to the ideal 
case, in particular plasma oscillations and surface plasmons 
(sometimes called surface plasmon polaritons). These are col- 
lective electron density waves with energies huj^ around ten 



electron volts (at typical metallic densities). These waves 
can be quantized and since huOp is larger than any exper- 
imentally relevant thermal energy, one can safely consider 
that bulk plasma modes are in the ground state |19|. This 
is not quite true for the surface plasmon modes that are con- 
fined to the surface of a metallic mirror. Their electronic 
excitation is accompanied by an electromagnetic field mode 
that is evanescent inside the cavity 1201 . Surface plasmons 
play an important role in many fields of physics. Let us 
only mention the plasmon-assisted light transmission through 
metallic structures ||2T]|22||23l, or dispersion forces between 
electronic Wigner crystals that are relevant for biomolecu- 
lar physics |24|. More generally, evanescent electromagnetic 
waves have a strong impact on the Casimir-Polder interaction 
between an atom and a surface as well as on the interaction 
between two surfaces at differences temperatures 1251 [261 . 

It is well known, indeed, that the Casimir effect, at short 
distances, is dominated by the coupling between the surface 
plasmons that propagate on two metallic mirrors. This has 
been pointed out in 1968 by Van Kampen and co-workers 1271 
who computed the Casimir energy for L <C Ap in terms of 
quasi-electrostatic (or non-retarded) field modes. In this limit 
the Casimir energy becomes ir7l [T2l[28l 
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-E^cas with a 1.790. 
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which is smaller than Eq.Q. Observe the different power law 
and the non-universal behavior as the result depends on the 
material parameter Ap. For metals used in modem experi- 
ments, Ap lies in the sub-micron range (107nm for Al and 
137nm for Cu and Au). This short-distance regime has been 
studied in much detail since Van Kampen 's paper, investigat- 
ing, for example, materials with a nonlocal response |29, 30l. 

As the mirror separation increases, retardation has to be 
taken into account, and Van Kampen 's result calls for a gen- 
eralization. This has been done by Schram in 1973 |31|, 
improving on a previous paper by Gerlach 1321 . Schram 
considered mirrors described by a non-dissipative dielectric 
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function, found the electromagnetic modes vibrating between 
these mirrors, and got the Casimir energy by summing their 
zero point energies. Among these modes, we find the retarded 
version of van Kampen's surface plasmon modes. Schram did 
not analyze separately their contribution and focused on the 
total energy, using a calculation based on the argument prin- 
ciple. Summerside and Mahanty investigated the joint effect 
of retardation and nonlocality on the surface plasmon modes 
at short distances |30|. 

In this paper we investigate more closely the influence of 
surface plasmon modes on the Casimir energy, covering both 
the non-retarded and retarded domains. This permits to ex- 
plore the experimentally relevant distance range around one 
micron where current precision experiments are performed. 
The plasmon modes are identified in a natural way in the sum 
over electromagnetic modes of the real cavity. We have shown 
previously that they have peculiar properties 1 33 1: one of them 
is purely evanescent while the dispersion relation of the other 
one changes its character from evanescent to propagating in- 
side the cavity (it crosses the light cone). In addition, the com- 
bined plasmonic contribution to the Casimir energy has the 
peculiarity to change sign as a function of distance L. Here, 
we derive and expand on these results in more detail and ex- 
hibit closed- form expressions valid at all distances. The main 
idea is to perform a re-parametrization of the dispersion re- 
lations that permit to evaluate analytically the relevant inte- 



grals. We recover van Kampen's result at short distances and 
discuss explicitly the asymptotic behaviour in the long dis- 
tance domain where retardation plays an important role. This 
regime was not covered in a previous paper by one of us ll34ll 
that performs an analysis of surface plasmons in the short- 
distance (non-retarded) regime. The analysis of the "photonic 
modes" (corresponding to waves that propagate in the cavity) 
will be the object of a following paper. For simplicity, we 
restrict here to zero temperature, the generalization to finite 
temperature being straightforward. 

The material is organized as follows. The basic method 
and the cavity modes are introduced in Sec. II. The disper- 
sion relation of the plasmonic modes is analyzed in Sec. Ill 
and Appendix [A| and their contribution to the Casimir energy 
given in Sec. Ill A. The Sees. Ill B and C discuss the short and 
large distance regimes. Our analysis concludes with a discus- 
sion of the sign of the Casimir interaction (Sec. HID) and of 
alternative splittings of the plasmonic dispersion relations that 
appeared recently in the literature (Sec. HIE). 



II. CASIMIR INTERACTION AND REAL CAVITY MODES 

In 1973 Schram proved the following mathematical identity 
L31J , exploiting the argument principle [5J 
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The left-hand side has the same structure as Casimir' s sum 
over zero point energies, but in this case the relevant modes 
are those of the real cavity. The notation [• • • ]^^^ signifies 
the difference of the expression in brackets for finite and in- 
finite mirror distance L. The right-hand side is nothing but 
the Lifshitz formula for the Casimir energy. Let us recall that 
Lifshitz adopted in 1955 |7| a fairly different viewpoint and 
computed the force as the average of the Maxwell stress tensor 
inside the cavity. He considered the electromagnetic fields as 
being radiated by fluctuating sources in the medium compos- 
ing the mirrors, similar to London's derivation of the Van der 
Waals force between atoms and molecules. The main point 
of Ref. |31 1 was to show that the Lifshitz approach yields the 
same result as the Casimir sum over zero-point energies, pro- 
vided the mirrors are non-dissipative. This is the case we fo- 
cus on here. 

The modes in Eq.([3]) are labelled by their polarization /i = 
TE, TM and the wavevector k = (kx^ky) parallel to the mir- 
rors; the perpendicular wavevector kz is defined in Eq.([6]) be- 
low. The are the reflection amplitudes that we take the 



same for both mirrors. The mode frequencies cj^(k) are re- 
lated to the zeros and the branch cuts of | 31 1 

D^[u;;k] = l-ri^[Lofe^"'^'^. (4) 

We adopt here the Fresnel formulas for the reflection ampli- 
tudes that for the case of thick mirrors read 1 36 1 

TE ^ ~ TM ~ e[(U;jAi: 

r — ^ — , r — r~r~r 

where 

kz = m = i\J |k|^ — uo'^ /c^ (6a) 

= ^J\kf -e[uj]coyc^ = ^K^+ul/c^ (6b) 

We choose signs for the square roots such that Re > 
and Im [hZi] < in Im [uj] > 0. This analytical continuation 
entails that Eq.Q has no solutions in the upper half plane. 
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Finally, e[cj] is the dielectric function; in the case of a metal 
the simplest description is given by the plasma model 



1.2 



eM = 1 



(7) 



where ujp is the plasma frequency, a constant which can be 
related to the specific physical properties of the metal. Up 
to uj ^ ujp the dielectric constant differs from unity so that 
the metal behaves different than the surrounding vacuum. For 
CO cOp the dielectric constant approaches unity and the metal 
becomes transparent. This is the way the plasma model imple- 
ments the high-frequency cutoff for the mirror reflectivity. 
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FIG. 1: Dispersion relations for TE-polarized modes between two 
metallic mirrors described by the plasma model (solid line), com- 
pared perfect conductors (dashed line). Mode frequency a;(k) and 
wavevector in the mirror plane, |k|, are normalized to the plasma 
frequency oOp. Mirror distance L = 1.75 Ap. The (blue) diagonal 
line is the light cone below which the field is evanescent in the cavity 
(evanescent modes). Above the thick solid line, the field propagates 
through in mirror material (bulk modes). 

In this model we neglect all the dissipation phenomena and 
we impose a local response to the electromagnetic field 1361 . 
From a physical point of view, it is a poor approximation to 
real metals at low frequencies (dissipation and non-locality, 
i.e., the anomalous skin effect are predominant) and high fre- 
quencies (absorption from intraband transitions). But at any 
rate, its mathematical simplicity allows explicit calculations 
to be pushed very far and to understand important physical 
behaviors. We are going to see that our principal result corre- 
spond to a frequency range high enough for the plasma model 
to be a good description of the metal. Let us stress, however, 
that the choice of the plasma model is the strongest approxi- 
mation we make and, within this model, all results we discuss 
are exact. 

The introduction of the dielectric properties of the mir- 
rors leads to a series of important modifications for the field 
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FIG. 2: Same as Fig[T]for TM-polarized modes. The solid red curves 
represent the plasmonic modes and the black curves the photonic 
modes. Note that one of the plasmonic modes crosses the light cone. 



modes. First of all, even in the simplest case, the plasma 
model, the dispersion relations cj^(k) cannot be written in 
terms of elementary functions. The results of a numerical cal- 
culation are shown in figs. [T]and|2](see details below). As we 
can see, imperfect reflection modifies the dispersion relation 
(solid lines) compared to a perfect reflector (dashed lines). 
We can distinguish three regions starting from above: 
Bulk modes occur for uj > ujbO^) = (^p + c^|kp)^/^ 
(shaded above the thick line); they propagate both in the cav- 
ity and inside the mirrors. These modes form a continuum that 
is mathematically represented by a branch cut of Eq.(|4]) in the 
complex cj-plane. This has to be taken into account carefully 
when applying the argument theorem |37|. The associated 
difficulties have led Schram to work instead with a mirror of 
finite thickness d where the continuum discretizes 1 31 J . For 
simplicity, we take here the limit of thick mirrors. 

Propagating (ordinary) cavity modes: they occur in the 
region above the light cone and below the bulk continuum, 
c|k| < uj < cjB(k). These modes are guided between the 
mirrors (note that the latter behave like a medium optically 
thinner than vacuum, < e[uj] < 1), leading to a discrete set 
of mode frequencies for a given k. In this region, the reflection 
coefficients ^ have unit modulus and a frequency-dependent 
phase. This leads to a shift of the cavity modes relative to 
perfectly reflecting mirrors, as is visible in Figs Tp 



Evanescent modes lie below the light cone, u < c|k| 
(shaded below the diagonal), and are the main focus of this pa- 
per. Their electromagnetic field exponentially decreases when 
going away from the vacuum-mirror interface, while it is al- 
lowed to propagate along the interface. Evanescent fields are 
of great interest in near field optics because they provide the 
link to sub- wavelength topographic features of a surface. In 
the context of the Casimir interaction, they are often underes- 
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timated due to their damped nature. We show here, however, 
that their contribution is all but a small correction, even at 
large distances 1331 . From a mathematical point of view, the 
optical properties of evanescent modes (reflection and trans- 
mission amplitudes) can be obtained from ordinary modes by 
a well-defined analytical continuation procedure 1351 . Solving 
the dispersion equation ^ in the evanescent sector, one finds 
two nondegenerate mode frequencies in only one polarization, 
at least for non-magnetic media. These modes are called "sur- 
face plasmons" (or "surface plasmon polaritons") l2Qll38ll39l . 
Their field amplitude decays exponentially away from the in- 
terface, and is associated with oscillating surface charge and 
surface current densities, as required by the equation of con- 



tinuity (see Fig|3]). On an isolated interface, surface plasmons 
correspond to the pole of r [u;] ; they occur when e[(jo] < — 1 
(i.e., CO < 00^/^2). For two interfaces, two surface plasmons 
exist and are coupled via their evanescent tails in the cavity. 
The resulting modes are given by the zeros of Eq.(|4]) for real 
hz and hZrn, and will be analyzed in detail in the following. 

Summarizing, we can see that for the TE-polarization, all 
modes lie above the light cone, while for TM-polarization, 
two modes enter the evanescent region in at least some range 
of wavevectors. We refer to these modes as "plasmonic"; they 
are the retarded generalization of van Kampen's coupled sur- 
face plasmon modes. Finally, we can re- write the Casimir 
energy as 



cavity modes 



bulk modes 



L^oo 



plasmonic modes (Epi) 




(8) 



photonic modes (-Eph) 



These contributions have no physical meaning on their own, 
i.e. one cannot measure them separately. The only observ- 
able is the total Casimir energy, which is the sum of all terms. 
However, evaluating them separately reveals striking features 
which suggest new possibilities to taylor the strength and the 
sign of the Casimir force. In the rest of this paper, we are go- 
ing to focus our attention on the plasmonic contribution Ep\ 
and shall discuss the remaining contributions to the Casimir 
energy in another paper. 



III. PLASMONIC MODES 



We plot again in Fig |4] the dispersion relation of the two 
modes in the first sum of Eq.([8]). They end up for large |k| 
below the light cone, i.e., the associated field is evanescent 
both in vacuum and in the mirrors. One branch that we call 
cc;_(k) lies entirely below the light cone. The second one, 
cj+(k) moves continuously into the cavity mode sector as k is 
decreased. The inset illustrates the smooth change in the spa- 
tial mode function. This mixed character justifies the name 
"plasmonic" that we use for both modes in the following. We 
discuss in Appendix [A| some general features of their disper- 
sion relations that can be obtained explicitly despite the fact 
that we have to deal with implicit functions. 



A. Contribution to the Casimir energy 

The plasmonic contribution is defined as the first sum on 
the r.h.s. of eq.([8]), namely 



pi 



E 



(9) 



Both modes tend to ujq (K) for L ^ oo so that we subtract the 
zero-point energy for two isolated surface plasmons. We are 
thus measuring the interaction energy arising from the cou- 
pling between the surface plasmons. 

Replacing the k-summation by an integral and using the 



scaled variables introduced in ( Al ), we get 



firA C KdK 



(10) 

To check the convergence at large K, we use the parametriza- 



tion of Eqs.( A3 ) and find the estimate 



1 



-2K\ 



(11) 



provided K max(l, ft^). For further details, see the dis- 



cussion around Eq.(^8]). The difficulty in Eq.([TO]) is that the 
dispersion relations Q±{K) are only known implicitly in the 
general case. We now show that using the parametrization of 
Appendix |A| the integrand can be brought into an explicit and 
elementary form. 
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FIG. 3: A surface plasmon mode is associate with an electronic 
charge oscillation bound to the surface of a body. For a single body 
the associated electric field is evanescent and, for a plane interface, 
the plasmon can be excited only by approaching from the vacuum 
side a medium with higher index of refraction and illuminating the 
latter in total internal reflection. Approaching two surfaces, the two 
respective surface plasmons couple through their evanescent field 
tails. A frequency splitting occurs giving rise to two new modes, the 
plasmonic modes. The antisymmetric and the symmetric (co-) 
mode have higher resp. lower energy than the isolated (non-coupled) 
mode ooq. The Casimir force associated with a;+ is then an anti- 
binding force (repulsive) while the oo- modes contribute an attrac- 
tive force. The plasmonic Casimir force arises from the (distance- 
dependent) balance of the two contributions. 




FIG. 4: A plot of the plasmonic dispersion relations a;+(k), cJo(k), 
a;-(k), as function of |k| for L = 0.2 Ap (Ap = 27rc/cjp). Fre- 
quencies and wavevectors are scaled to the plasma frequency cOp and 
cjp/c, respectively. Inset: magnetic field amplitude for chosen points 
along the branch (k), as labelled by the circles. 



It is useful to scale the energy Eq.([TO| to the perfect-mirror 
Casimir energy E^cas [Eq.Q] 



(12) 



POO 

^Pl = -^ / E Cana{K)KdK (13) 

a=±,0 



for the plasmonic Casimir energy; note that it depends on the 
distance only via the dimensionless parameter ftp. 

For each of the branches ft a {K), we now change to the inte- 
gration variable z = {nL)'^ . The Jacobian (the prime denotes 
the derivative) 



dK'^ = 2KdK = dz^ 9a{z)dz 



with ga{z) defined in Appendix [a] leads to 

^ a=±,0 ^ 



(14) 



(15) 



The integration paths are now r+ = —z^ oo, and r_^o = 

. . . oo where z+ is defined in Eq.( [A8[ ). 

One of the two terms under the integral can be integrated 
immediately, leading to 



o=±,0 



Y \ ga{z)^fgjz)dz = ^ 



o=±,0 



^3/2 



(^) 



(16) 

where the function in brackets has to be evaluated at the end 
points of the respective integration domains. The upper limit 
contributions {z = oo) cancel under the subtractions. At the 
lower limit, g^^Q{0) vanishes because the dispersion relations 
reach u = (see SeciAb. We are thus left with g^'^{—Z-^) = 

3/2 
Z^ . 

Putting the propagating sector of the mode ft-^{K) into 
a separate integral, the correction factor for the plasmonic 
Casimir energy can be rewritten as 



/ '^Ca\/ga{z)dz^ / ^/g^ 

"^0 a '^-^+ 



{z)dz 



2 3/2 

(17) 



In the first integral, the functions ga{z) are real. For 2: ^ oo, 
the functions g±{z) approach go{z) exponentially fast. An 
expansion in leads to 

E (^aVgJ^ « -n^e-'^^f{z/Ql) (18) 



with c+ = C- = 1, Co = —2. We call T^pi the correction factor 



where the function f{z/ft'^) is bounded and tends to 1/ (4\/2) 
for z l]p. This secures the convergence at large z of 
the first integral in Eq.(p7|). The second integral is finite be- 
cause g-\-{z) is bounded and the integration domain is finite 
[see Eq.( |A9| )]. Both the second integral and the third term 
in Eq.(p^ are related to the propagating segment of the plas- 
monic mode ft-^{K). 

The great advantage of Eq.([T7]) compared to Eq.([TO]) is that 
now the integrands are expressed in terms of simple analytic 
functions and there is no need to integrate implicit functions 
whose evaluation is only possible numerically. We also gain 
for analytical calculations since the discussion of the distance 
dependence (via the parameter ftp oc L/X^) can be done in 
a transparent way. We show in the following that one gets 
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asymptotic expressions for small and large values of ftp, the 
only variable on which the correction factor r/pi depends after 
the integration. 

Fig|5] shows a plot of r^pi as function of L/Ap = Q^/ (27r). 
Note the increase linear in L for small distances and a sign 
change at large L, with a power law ex L^/^. In the next 
two sections we analyze these limits analytically. At short 
distance r]pi reproduces exactly the correction factor known 
for the total Casimir energy. 



B. Short distance asymptotics 

The distance enters the correction factor r]pi [Eq.(p7|)] via 
the dimensionless parameter l^p, and we get the short-distance 
asymptotics in the limit fip ^ 1. This has been discussed in 
previous papers ||28l|32j[34||40l, but the asymptotics turns out 
to be tricky at next-to-leading order. 

The first order expansion in of the functions ga{z) 
yields EH 



(19a) 



0.08 



Qp L 

where the numerical constant a ^ 1.790 arises from 



a 



"71 



\/l + e-v^ + Vl- e-^ - 2^ dz (19b) 



The separate contributions of the modes fl-^{K) and fl_{K) 
are a+ ^ —12.225 (repulsive) and a_ ^ 14.015 (attractive). 
The plasmonic Casimir energy in this regime thus scales like 
Ahujp/L'^ and is reduced compared to the perfect mirror case 
(V « 1) 141 1. 

The contribution of the propagating part of cc;+ is of the 
third order in [see Eq.(|A9|)] 



g+[z)dz > 9+{^)z+ 



^+0^p 



(20) 



and can therefore be neglected. The same argument holds for 
the term (z+)^/^ l^p. In other words, the plasmonic contri- 
bution comes essentially from the evanescent sector {z > 0). 



We note that the result ( |19a| ) yields exactly the short-distance 
behavior of the full Casimir energy which is thus dominated 
at short distance by the interactions between surface plas- 
mons Il28l 132113411401. 

Note that Eq.( |19a| ) follows from an expansion of the ga{z) 
to first order in ft^. It is worth stressing that this expansion 
scheme does not work at higher orders, the series being an 
asymptotic one and not uniformly convergent. Each integral 
obtained by this method at higher orders is divergent, except 
the first one given in Eq.( |19a| ). To avoid this problem, we use 
an alternative method and write the functions ^± (^) as follows 



g±{z) = go{z) 



(21) 




0.02 0.04 0.06 0.08 
L/Ap 

FIG. 5: A plot of the plasmonic Casimir energy, normalized to the 
perfect mirror case (i.e., the plasmonic correction factor ?7pi) vs. the 
scaled distance L/Xp. Light gray solid line: approximation ( [24| . 
The inset illustrates the short-distance behaviour beyond the linear 
term: we plot [?7pi — a{L/Xp)]/Qp (black line) and compare to the 
expression a + 61og(27rL/Ap) in Eq.([24]) (light gray line). 



with 



P 



_^^ go{z)-nl/2 



(22) 



Now for z > 0, IpI is bounded by unity and decays rapidly to 
zero as z ^ oo [i28l . To compute the integral of [g±{z)]^^'^, 
we expand in powers of p and get the series 



[9±iz)] 



1/2 _ 



5o(z)(l ± e- 



{±pr 



^r(i-n) n! 

n=0 ^ ^ ^ 



(23) 

where r(-) is the gamma function. Taking the n = term, 
the integration over z leads to ( |19a| ). Higher order terms can 
be calculated explicitly, but the resulting expressions are cum- 
bersome and will not be reported here. Including the next-to- 
leading order terms, we find 



(24) 



where a ^ 0.63 and b = ^ 1.026. This is plotted as 

gray line(s) in Fig|5] the inset providing a zoom on the cubic 
and logarithmic terms (see caption). Note that the term aQ^ 
gives a distance-independent correction to the Casimir energy 
and cancels when the force is computed. The presence of the 
logarithmic correction is due to the non-uniform convergence 
of the asymptotic series. We find from ([24]) that the Casimir 
force does not feature a logarithmic correction at short dis- 
tance. 
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C. Large distance asymptotics 

The curve r]pi{L) in Figjs] shows that the plasmonic mode 
contribution is negative (repulsive) at distances L > 0.08 Ap. 
Mathematically, this can easily be seen from the large ftp 
asymptotics of r^pi. One can check that the integrand of the 
first integral in Eq.([T7]) is significantly different from zero only 
for z ~ 1. This suggests the following expansion of the ga{z) 
for l^p > 1 



(25a) 

(25b) 
(25c) 



Moreover, the expansion to leading order in ft'^ \z\ can 
also be performed in the integral over the propagating sector 
in ^Tl) because the integration domain is limited to — tt^ ^ 
Z-^ < z < 0. Finally, we find that the integrated term in 
Eq.(p^ gives a negligible contribution so that to leading order. 



L40J with 



,3/2 



coth[|] + Wtanh[|]-2)d7/ 



(26) 



This expression can be evaluated numerically, giving as result 
r = 29.75 (i.e. the sum of 8.90 (+ mode, evanescent sec- 
tor), — 7.23 (— mode, evanescent sector), and 28.09 (+ mode, 
propagating sector). Note the large contribution of the prop- 
agating segment and the near cancellation of the two evanes- 
cent branches. 

Since r^pi is negative at large distances, the plasmonic con- 
tribution provides a repulsive contribution to the Casimir in- 
teraction that scales like -\-Ahy/uJ^/L^^'^. This is balanced 
in the total Casimir energy by the contributions of photonic 
modes (cavity and bulk modes), recovering the attractive 
large-distance power law ^cas —Ahc/L^. 



D. Cancellations and signs 

We conclude our analysis by suggesting an interpretation of 
the signs of the plasmonic contributions to the Casimir energy. 
It is clear that Ep\ is due to the shift in the plasmon mode 
frequency relative to the isolated interface [see Eq.([T0])]. 

This can be also interpreted as a reshuffling of the density 
of modes due to the coupling by the interface, the total number 
of modes remaining constant. To make this more quantitative, 
we re-write the plasmonic Casimir energy as 



k,a=± 



^ J duj uj [pa{u;) - po{uj)] 
(27) 



a=± 



where the mode densities are defined as usual by (a = ±, 0) 

Pa{u;) = Y,^i^-^am^ (28) 

k 

that depend on the distance L for a = ±. Note that the u in- 



tegral in ( 27 ) does not converge if taken over the pa (uj) alone. 



This is due to the flat large- /c asymptote of the plasmonic dis- 
persion relations. More explicitly, the density of modes can 
be calculated as 



Aka{u;;L) 
27r 



dka{u;;L) 



(29) 



where, ka{uj) is the inverse function to uJaO^), and the deriva- 
tive is just the inverse group velocity at a given frequency 
UJ. We find a behaviour Pa{^) oc (cjgp — cc;^)~^ when uj 
approaches the asymptotic value uj^p = ujp/V2 of the dis- 
persion relation (the surface plasmon resonance in the quasi- 
static limit). This peak is exactly cancelled in the difference 
^P±{^) = P±{^) — Po(^) that we plot in Fig|6]for a given 
distance L. The precise behaviour of the curves changes with 
the distance (at smaller L, for example, P-\-{uj) is nonzero 
for UJ > ujsp), but the following qualitative features are sta- 
ble, (i) The mode cj+(/c) shows a gap between and t^+(0), 
and the difference dp-^{uj) is only due, for uj < uj-^{0), to 
the subtracted isolated surface plasmon (dashed line). Just 
at this frequency, the mode density P-\-{uj) jumps to a posi- 
tive value. This behaviour is due to the quadratic shape of the 



lower band edge in uj-^{k). As uj 



3, 6p-^{uj) > because 



uj-^{k) is shifted upwards relative to uJo{k) (the group veloc- 
ity is smaller). This mode is hence an 'anti-binding one' 1341 . 
(ii) The mode uj-{k) has a linear dispersion for small k, and 
the difference in mode density can be worked out as the pos- 
itive quantity Sp-{uj) ex uj/Qp ex ujL (dashed line). This 
mode is hence anti-binding in this region as well, (iii) Near 
the frequency uj^p, the mode uj-^{k) {k)] gives a repulsive 
[attractive] contribution to the Casimir energy, respectively. 
Summing over both modes yields a repulsive or attractive re- 
sult depending on L, because the relative weight of the bind- 
ing and anti-binding regions changes. Coming back to fre- 
quency shifts, it is easy to see from the large K expansion of 
Eqs.(|A3|) that the following inequality holds 



|cj_(k,L) -cjo(k)| > |cj+(k,L) -cc;o(k)| (k>cc;p/c), 

(30) 

At short distance, the plasmonic Casimir energy (which is ac- 
tually the total Casimir energy) is thus attractive, as is well 
known (see also EqfTS]). 

Let us finally note that as one moves away from the large- 
k regime, retardation becomes increasingly important. The 
change in sign of the plasmonic Casimir energy can thus be 
seen as well as a consequence of the finite speed of light. 



E. Cutting the mode branch 

Recently, there has been some discussion on the way 
to split the field modes into photon-like and plasmon-like 
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parts BOl 1461 . We comment in this section on the numbers 
one can obtain when the plasmonic mode ooj^{k) is segmented 
in a different way. (The mode ij-{k) is subject to no con- 
troversy.) The main conclusion we draw from this discus- 
sion is that the large distance behavior is dominated by mode 
branches near the light line. In addition, the sign is sensitive to 
the chosen subtraction (renormalization), and it may happen 
that under this procedure, a pure evanescent branch ends up 
being counted among photonic modes. We also suggest that 
the branch of the plasmonic mode ooj^{k) that enters the prop- 
agating sector is perhaps one of the best examples of Casimir 
repulsion due to a standing wave mode. Consider the corre- 
sponding pressure: it is repulsive due to photons bouncing on 
the mirrors. The attractive force for a perfect cavity arises, all 
things told, from the subtraction of a similarly repulsive pres- 
sure from a standing wave mode continuum (reflected from 
the mirrors' backfaces). Now, the counterpart for the plas- 
monic mode is a single-interface evanescent mode with zero 
pressure so that the repulsive force survives the subtraction. 



Bordag BOl is calling 'plasmon mode' only the evanes- 
cent branch of ooj^(k) that exists for A: = |k| > = 
ujp/ {cy^l + (see Fig|7] top). The segment within the 

light cone actually does not appear explicitly in Eq.(24) of 
Ref. BOl , but is implicitly contained in the total Casimir en- 
ergy (the photonic contribution is computed by subtracting the 
plasmonic one). The evanescent segment of cc;+(/c) is renor- 
malized by subtracting the isolated surface plasmon, cjo(^), 
over the same range kc < k < oo, sls shown in Fig |7] (top). 
The range < k < kc is left out (although it depends on L 
via kc). This subtraction is sufficient to get a vanishing en- 
ergy as L ^ oo because uj-^{k) ujqO^) exponentially fast 
for k > kc. (In addition, kc 0.) The integration over the 
branches chosen in Ref.|40| corresponds to the following cor- 




FIG. 7: Illustration of different segmentation of the plasmonic modes 
and the chosen renormalization. Tfiick lines mark tfie segments tfiat 
are taken into account in tfie different approacfies. We write oopr and 
iOev for tfiose parts of tfie mode wfiere tfie field between tfie 

mirrors is propagating or evanescent, respectively. 
Top: Bordag |40|, tfie modes u;+,o(A:) (red, blue) start at tfie 
wavevector kc wfiere a;+ (k) reacfies tfie ligfit cone. Middle: one pos- 
sibility suggested by tfie comment of Lenac |46|. Tfie mode iOev{k) 
is continued, for < k < /cc, by tfie ligfit line uj — ck (red) 
and renormalized by tfie entire brancfi of cjo(^) (blue). A particu- 
lar splitting of tfie Lifsfiitz formula into propagating and evanescent 
modes turns out to yield tfie same result. Bottom: anotfier possi- 
bility compatible witfi Lenac 's paper. Only tfie evanescent brancfi 
^ew{k) (k > kc, red) is taken into account and renormalized by 
ooo{k) (k > 0, blue). 



rection factor to the Casimir energy: 

Zc 

~ 2 





y^dz+-{nl{K,)-K^,) 



(31) 



with 



^ ^ a=±,0 



(32) 
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Here, = kcL and Zc solves the equation = go{zc) 
(at this parameter value, the dispersion relation u;o{k) reaches 
k = kc). We have checked that at short distance, this correc- 
tion is negligible compared to the leading order r/pi oc cOpL. 
At large distance, ho wever the integrals in ( [3T] ) are both of or- 



waves. With the change of variable k ^ z = (/^L)^, we see 



der l^y^ (see Section llIIC|) and the difference l^o(^c) - ' ^^^^ ^^^^^ straightforwardly to Eq.([32 ) 



that the two factors in the second line of p4\ have simple ze- 
ros (poles) at the mode frequencies u;± (cjq), respectively [see 
Eq.( |A3a| )]. A calculation using the argument principle and the 
symmetry property e(— cj) = e{u;) for lossless response func- 



too. Their contributions come with different signs, leading in 
the end to a correction factor that is attractive and scales like 

1/2 

77b ^ 1.6240 at large distance. 

A similar analysis can be done for the mode definition 
sketched in Fig|7] (middle): the plasmonic mode is contin- 
ued along the light line for k < kc and renormalized by the 
entire dispersion branch uJo{k). For L ^ oo, as ^ 0, 
the renormalized energy vanishes. The corresponding cor- 



rection factor is given by t/l [Eq.(32)] which does not con- 
tain any integrated term. The short-distance behaviour is the 
same as in the present paper, and at large distance, we have 

1/2 

r]h{L) ~ —1.6600Qp . This corresponds to repulsion as with 
our convention, but with a smaller numerical coefficient. We 
argue below that this result can also be obtained by a splitting 
of the Lifshitz formula for the Casimir energy. Let us men- 
tion that if the segment < k < kc of the light line is not 
taken into account (FiglT] bottom)) then the large distance be- 

3 /2 

haviour shows an attractive term r] ex . Both results do 
not fit with the curves presented by Lenac (461, although our 
calculation tries to follow the spirit of his description. It is not 
clear to us from his sparse description which renormalization 
scheme was used in the end. 

The Lifshitz approach to the Casimir energy leads to the 
correction factor t^l as follows. We write the right-hand side 
of Eq.([3]) in the equivalent form 



-Im 



oc 



^nc.£logI^,[c.,k] 



(33) 



where the dispersion function I)^ [cj, k] is defined in This 
expression has a structure very similar to the so-called "argu- 
ment principle" where the zeros (and poles) of the argument 
of the logarithm define the eigenfrequencies of the system (of 
the reference system), respectively |5|, and each mode con- 
tributes its zero point energy. In other words, the imaginary 
part of the logarithmic derivative can be read as a density of 
modes (suitably renormalized). We isolate the contribution 
of evanescent modes by restricting the cc; -integration domain 
to < < c|k| (so that = ^|kp — co'^/c^ is real as it 
should for evanescent waves). As discussed in Sec|ll| zeros 
and poles of D^[uj^ k] occur for the plasma model only for 
TM-polarized evanescent waves. A simple calculation leads 
to 



1 - {rr're-'^'^ = (1 - {r^^fe-'^"^) x (34) 



where the functions g±{z) defined in ( A3b[|A3c| ) appear. The 
factor involving r shows no singularities for evanescent 



IV. CONCLUSION AND DISCUSSION 

In this paper we evaluate the contribution of plasmonic 
modes to the Casimir force using the plasma model to de- 
scribe the optical response of the medium. Simple analyt- 
ical expressions are found, in particular for the small and 
large distance asymptotics. We introduced a correction factor 
r]p\{L) that gives the plasmonic contribution to the Casimir 
energy, Epi{L), in units of the Casimir energy Ecas{L) oc 
— [Eq.M]. It turns out that r]pi{L) is small but pos- 
itive at short distance, correctly reproducing van Kampen's 
result IZTTI . Quite surprisingly, the plasmonic contribution 
changes sign at the fairly short distance L/Xp ~ 0.08. For 
larger cavity lengths, r]p\{L) becomes negative and leads to the 
unusual scaling Ep\ ex asL oo. This behaviour 

clearly shows that the plasmonic modes are much more im- 
portant for the Casimir effect than usually anticipated. They 
do not only dominate in the short distances limit, but they also 
give a large repulsive contribution at large distances. 

We have calculated as well (see also [[40|) the photonic 
mode contribution that turns out to be a monotonous function 
of the distance L (Fig [8]); it actually approaches a constant as 
L ^ 0. Its large distance behaviour contains as leading or- 
der a negative term that exactly cancels the plasmonic 
contribution. The Casimir energy is thus the balance of two 
contributions of equal magnitude which nearly cancel each 
other. 

It would be interesting to investigate if a change in the field- 
mirror coupling could somehow influence this detailed bal- 
ance and therefore the value or even the sign of the Casimir 
force. This could be the case for nano structured surfaces, 
since the plasmonic modes are associated with the electron 
charge density oscillations at the vacuum/metal interface. 
This route has already been explored within a different con- 
text, that of metamaterials in the visible frequency range. It 
has been shown that arrays of metallic dots or rods Il42ll43ll44ll 
exhibit a strong magnetic response in the visible, including 
a band with negative magnetic permeability. This behaviour 
arises again from plasmon modes: they are here concentrated 
on the metallic particles and their characteristics can be tuned 
with the particle shape. In the array, the plasmons delocalize 
and lead to a resonant electric and magnetic response. A sig- 
nificant modification of the Casimir force, even a change in 
sign, could be realistic with these materials | 45 1. 

The limits of validity of our results are imposed by the ap- 
plicability of the plasma model. The approach is not intended 
to make quantitative predictions because the response of intra- 
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FIG. 8: A plot of plasmonic ^pi, photonic and total Casimir E 
energy vs. distance L, normalized to the plasma wavelength Ap. We 



normalize the energy to (27r) 



hcTv-^A 
720A3 • 



The plasmonic energy shows a 



maximum for L ^ 0.16 Ap (the corresponding force changes sign), 
while the photonic energy provides an attractive contribution at all 
distances. 



band transitions in real metals would require a more compli- 
cated dielectric function. We are also restricted to fairly short 
distances where the relevant mode frequencies are sufficiently 
large compared to the dissipation rate. From FigjSj one can 
see, however, that the most striking effects due to plasmonic 
modes (change in sign and near cancellation of plasmonic 
and photonic contributions) indeed appear at short distances, 
L < Ap, where the lossless plasma response is a suitable ap- 
proximation. We therefore believe that at least for this range 
of distances our results are quite generally valid. 

In our description, the main responsible for a repulsive 
Casimir interaction is the plasmonic mode cj^. This mode 
crosses the border between the evanescent sector and the prop- 
agative sector, and we have it considered as being completely 
part of the 'plasmonic' set of modes |33|. This is an 'adia- 
batic' definition that is strongly suggested by the continuous 
change in the mode function plotted in Fig|4] Other splittings 
into evanescent and photonic modes have been applied in the 
literature |[40| |46l, and we have given a brief review in Sec- 
tion HIE The total Casimir energy is of course immune to 



these wordings. However, if one considers a structured sur- 
face, the mode branch cj+ will change as a whole, and by ana- 
lyzing this change, one could easily predict the corresponding 
modification of the Casimir energy. 
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APPENDIX A: PLASMONIC DISPERSION RELATIONS 

We collect in this appendix some properties of the solutions 
oo±(k) of the dispersion relation Eq.(|4]) in the evanescent sec- 
tor. Similar discussions can be found in Refs. ll38l[39l . 

To simplify the notation, we use the dimensionless vari- 
ables 



i^=|k|L, n = ujL/c, 



(Al) 



and write the ratio L/Ap = Qp/27r. By symmetry, the mode 
frequencies do not depend on the angle of k within the mirror 
plane. The zeros of Eq.(|4]) and can be split into 'odd' and 
'even' cases 



(A2) 



For the TM-polarization, the solutions to these equations can 
be written in the parametrized form 



where 



(A3a) 



(A3b) 



(A3c) 



The parameter z = {nL)'^ runs over a part of the real axis such 
that z^g±{z) > in Eq.( [A3a ): it is positive in the evanescent 
sector and negative in the propagating one. The square root is 
chosen with Im < 0. For completeness, we give here the 
form of with negative argument: 



-^/tan(^V^) 



(A3d) 



it takes positive values as long as z' is not too large (see 
Eqs.( |A8|A9| ) below). 

We shall also need the surface plasmon dispersion for a sin- 
gle metal- vacuum interface, VtQ{K). It can formally be found 
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by taking the limit L ^ oo (while hz > 0) in Eqs.( |A3b|A3c| ), 
leading to the function 



9o{z) 



(A3e) 



that parametrizes according to the scheme of Eqs.(A3a) the 
poles of [Vt] . We note for completeness the explicit ex- 
pression for the single interface 



^0 W = \ 



(A4) 



(This is of course independent of L since the scale factor 1/L 
can be removed throughout.) 

In the evanescent sector, both functions z + g±{z) are 
monotonous so that exactly one evanescent solution is found 
for a given K. The limit z ^ co corresponds to large wave 
vectors K where we have 



K- 



±,0. 



(A5) 



The two surface plasmons decouple in this limit. A more 
detailed analysis of the asymptotic behaviour is given in 
SecHnAl 

The modes Q-{K) and Qo{K) remain in the evanescent 
sector fovK 0: analyzing the limit z ^ 0, one gets the 
linear dispersions 



K 



K 



(A6) 
(A7) 



The mode Q.j^{K) crosses the light line when z = 0, 
corresponding to the frequency l^+o = [^+(0)]^^^ = 
(^p/y^I + l]p/2. The dispersion relation enters the cavity 
sector in a continuous and differentiable way. The same is 
true for the associated mode function, as illustrated in Fig|4] 
(inset). It reaches i^^ = at a parameter value z = —Zj^ that 
is the solution of z + g+{z) = 0. This can be simplified to the 
implicit equation 



= l]p cos(^y^) 

from which one gets the exact inequalities < 
min(l]p, TT^) and the asymptotics 



for ^Ip <C TT 
for TT 



(A8) 



< 




(A9) 



In the first case, the mode ends at the bulk mode continuum 
(Fig|4]is close to this situation). The second case implies that 
at large distance (L Ap), the mode stays close to the light 
line, as in Fig|2] One can also calculate that the generic be- 
haviour of the dispersion relation for small K is quadratic: 



(AlO) 
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